THE EFFECTS OF PARAMETRIC UNCERTAINTIES IN SIMULATIONS OF A REACTIVE PLUME USING A LAGRANGIAN STOCHASTIC MODEL

A combined Lagrangian stochastic model with micro mixing and chemical sub-models is used to investigate a reactive plume of nitrogen oxides (NOx) released into a turbulent grid flow doped with ozone (O3). Sensitivities to the model input parameters are explored for high NOx model scenarios. A wind tunnel experiment is used to provide the simulation conditions for the first case study where photolysis reactions are not included and the main uncertainties occur in the parameters defining the turbulence scales, the source size and the reaction rate of NO (nitric oxide) with O3. Using nominal values of the parameters from previous studies, the model gives a good representation of the radial profile of the conserved scalar [NOx] compared to the experiments, although the width of the simulated profile is slightly smaller, especially at longer distances from the source. For this scenario, the Lagrangian velocity structure function coefficient has the largest impact on simulated [NOx] profiles. At the next stage photolysis reactions are included in a chemical scheme consisting of eight reactions between species NO, O, O3 and NO2. The high dimensional model representation (HMDR) method is used to investigate the effects of uncertainties in the various model inputs resulting from the parameterisation of important physical and chemical processes in the reactive plume model, on the simulation of primary and secondary chemical species concentrations. Both independent and interactive effects of the parameters are studied. In total 22 parameters are assumed to be uncertain, among them the turbulence parameters, temperature dependant rate parameters, photolysis rates, temperature, fraction of NO in total NOx at the source and background concentration of O3. Only uncertainties in the mixing time scale coefficient and the structure function coefficient are responsible for the variance in the [NOx] radial profile. On the other hand, the variance in the [O3] profile is caused by parameters describing both physical and chemical processes

THE EFFECTS OF PARAMETRIC UNCERTAINTIES IN SIMULATIONS OF A REACTIVE PLUME USING A LAGRANGIAN STOCHASTIC MODEL

A combined Lagrangian stochastic model with micro mixing and chemical sub-models is used to investigate a reactive plume of nitrogen oxides (NOx) released into a turbulent grid flow doped with ozone (O3). Sensitivities to the model input parameters are explored for high NOx model scenarios. A wind tunnel experiment is used to provide the simulation conditions for the first case study where photolysis reactions are not included and the main uncertainties occur in the parameters defining the turbulence scales, the source size and the reaction rate of NO (nitric oxide) with O3. Using nominal values of the parameters from previous studies, the model gives a good representation of the radial profile of the conserved scalar [NOx] compared to the experiments, although the width of the simulated profile is slightly smaller, especially at longer distances from the source. For this scenario, the Lagrangian velocity structure function coefficient has the largest impact on simulated [NOx] profiles. At the next stage photolysis reactions are included in a chemical scheme consisting of eight reactions between species NO, O, O3 and NO2. The high dimensional model representation (HMDR) method is used to investigate the effects of uncertainties in the various model inputs resulting from the parameterisation of important physical and chemical processes in the reactive plume model, on the simulation of primary and secondary chemical species concentrations. Both independent and interactive effects of the parameters are studied. In total 22 parameters are assumed to be uncertain, among them the turbulence parameters, temperature dependant rate parameters, photolysis rates, temperature, fraction of NO in total NOx at the source and background concentration of O3. Only uncertainties in the mixing time scale coefficient and the structure function coefficient are responsible for the variance in the [NOx] radial profile. On the other hand, the variance in the [O3] profile is caused by parameters describing both physical and chemical processes

Irreversible loss in marine ecosystem habitability after a temperature overshoot

Anthropogenic warming of the oceans and associated deoxygenation are altering marine ecosystems. Current knowledge suggests these changes may be reversible on a centennial timescale at the ocean surface but irreversible at deeper depths even if global warming were to ameliorate. In contrast, the marine ecosystem’s response to these persistent changes remains poorly elucidated. Here we explore to what extent global warming may drive alterations in marine habitats by exploring the evolution of a metabolic index that captures marine organisms’ ecophysiological response to both temperature and oxygen changes, throughout an idealised ramp-up/ramp-down atmospheric carbon dioxide concentration and an overshoot scenarios. Using a multi-model approach; we find that changes in ocean temperature and oxygen drive a centuries-long irreversible loss in the habitable volume of the upper 1000 m of the world ocean. These results suggest that the combined effect of warming and deoxygenation will have profound and long-lasting impacts on the viability of marine ecosystems, well after global temperatures have peaked.publishedVersio

RECCAP2 Future Component: Consistency and Potential for Regional Assessment to Constrain Global Projections

This is the final version. Available from Wiley via the DOI in this record. Data Availability Statement: All CMIP6 model output datasets analyzed during this study are available online at https://esgf-node.llnl.gov/search/cmip6/ and code required to reproduce figures is available at https://github.com/ChrisJones-MOHC/RECCAP2Future_2023 (ChrisJones-MOHC, 2023) and Zenodo at https://doi.org/10.5281/zenodo.8420250.Projections of future carbon sinks and stocks are important because they show how the world's ecosystems will respond to elevated CO2 and changes in climate. Moreover, they are crucial to inform policy decisions around emissions reductions to stay within the global warming levels identified by the Paris Agreement. However, Earth System Models from the 6th Coupled Model Intercomparison Project (CMIP6) show substantial spread in future projections—especially of the terrestrial carbon cycle, leading to a large uncertainty in our knowledge of any remaining carbon budget (RCB). Here we evaluate the global terrestrial carbon cycle projections on a region‐by‐region basis and compare the global models with regional assessments made by the REgional Carbon Cycle Assessment and Processes, Phase 2 activity. Results show that for each region, the CMIP6 multi‐model mean is generally consistent with the regional assessment, but substantial cross‐model spread exists. Nonetheless, all models perform well in some regions and no region is without some well performing models. This gives confidence that the CMIP6 models can be used to look at future changes in carbon stocks on a regional basis with appropriate model assessment and benchmarking. We find that most regions of the world remain cumulative net sources of CO2 between now and 2100 when considering the balance of fossil‐fuels and natural sinks, even under aggressive mitigation scenarios. This paper identifies strengths and weaknesses for each model in terms of its performance over a particular region including how process representation might impact those results and sets the agenda for applying stricter constraints at regional scales to reduce the uncertainty in global projections.European Union’s Horizon 2020European Union’s Horizon 2020European Union’s Horizon 2020Joint UK BEIS/Defra Met Office Hadley Centre Climate ProgrammeCarbonWatch-NZ Endeavour Research ProgrammeSão Paulo Research FoundationSão Paulo Research FoundationSão Paulo Research FoundationNational Science FoundationAndrew Carnegie Fellow ProgramCNPqKorea Ministry of EnvironmentNatural Environment Research Council (NERC)Natural Environment Research Council (NERC)National Environmental Science Progra

Is there warming in the pipeline? A multi-model analysis of the Zero Emissions Commitment from CO₂

The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory

Attribution of multi-annual to decadal changes in the climate system: The Large Ensemble Single Forcing Model Intercomparison Project (LESFMIP)

Multi-annual to decadal changes in climate are accompanied by changes in extreme events that cause major impacts on society and severe challenges for adaptation. Early warnings of such changes are now potentially possible through operational decadal predictions. However, improved understanding of the causes of regional changes in climate on these timescales is needed both to attribute recent events and to gain further confidence in forecasts. Here we document the Large Ensemble Single Forcing Model Intercomparison Project that will address this need through coordinated model experiments enabling the impacts of different external drivers to be isolated. We highlight the need to account for model errors and propose an attribution approach that exploits differences between models to diagnose the real-world situation and overcomes potential errors in atmospheric circulation changes. The experiments and analysis proposed here will provide substantial improvements to our ability to understand near-term changes in climate and will support the World Climate Research Program Lighthouse Activity on Explaining and Predicting Earth System Change.publishedVersio

Carbon-concentration and carbon-climate feedbacks in CMIP6 models, and their comparison to CMIP5 models

Abstract. Results from the fully-, biogeochemically-, and radiatively-coupled simulations in which CO2 increases at a rate of 1 % per year (1pctCO2) from its pre-industrial value are analyzed to quantify the magnitude of two feedback parameters which characterize the coupled carbon-climate system. These feedback parameters quantify the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate. The results are based on eight comprehensive Earth system models from the fifth Coupled Model Intercomparison Project (CMIP5) and eleven models from the sixth CMIP (CMIP6). The comparison of model results from two CMIP phases shows that, for both land and ocean, the model mean values of the feedback parameters and their multi-model spread has not changed significantly across the two CMIP phases. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The sensitivity of feedback parameters to the three different ways in which they may be calculated is shown and, consistent with existing studies, the most relevant definition is that calculated using results from the fully- and biogeochemically-coupled configurations. Based on these two simulations simplified expressions for the feedback parameters are obtained when the small temperature change in the biogeochemically-coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters allows identification of the reasons for differing responses among ocean and land carbon cycle models. </jats:p

Global Surface Ocean Acidification Indicators From 1750 to 2100

Accurately predicting future ocean acidification (OA) conditions is crucial for advancing OA research at regional and global scales, and guiding society's mitigation and adaptation efforts. This study presents a new model-data fusion product covering 10 global surface OA indicators based on 14 Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6), along with three recent observational ocean carbon data products. The indicators include fugacity of carbon dioxide, pH on total scale, total hydrogen ion content, free hydrogen ion content, carbonate ion content, aragonite saturation state, calcite saturation state, Revelle Factor, total dissolved inorganic carbon content, and total alkalinity content. The evolution of these OA indicators is presented on a global surface ocean 1° × 1° grid as decadal averages every 10 years from preindustrial conditions (1750), through historical conditions (1850–2010), and to five future Shared Socioeconomic Pathways (2020–2100): SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. These OA trajectories represent an improvement over previous OA data products with respect to data quantity, spatial and temporal coverage, diversity of the underlying data and model simulations, and the provided SSPs. The generated data product offers a state-of-the-art research and management tool for the 21st century under the combined stressors of global climate change and ocean acidification. The gridded data product is available in NetCDF at the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Information: https://www.ncei.noaa.gov/data/oceans/ncei/ocads/metadata/0259391.html, and global maps of these indicators are available in jpeg at: https://www.ncei.noaa.gov/access/ocean-carbon-acidification-data-system/synthesis/surface-oa-indicators.html

Carbon-concentration and carbon-climate feedbacks in CMIP6 models, and their comparison to CMIP5 models

Abstract. Results from the fully-, biogeochemically-, and radiatively-coupled simulations in which CO2 increases at a rate of 1 % per year (1pctCO2) from its pre-industrial value are analyzed to quantify the magnitude of two feedback parameters which characterize the coupled carbon-climate system. These feedback parameters quantify the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate. The results are based on eight comprehensive Earth system models from the fifth Coupled Model Intercomparison Project (CMIP5) and eleven models from the sixth CMIP (CMIP6). The comparison of model results from two CMIP phases shows that, for both land and ocean, the model mean values of the feedback parameters and their multi-model spread has not changed significantly across the two CMIP phases. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The sensitivity of feedback parameters to the three different ways in which they may be calculated is shown and, consistent with existing studies, the most relevant definition is that calculated using results from the fully- and biogeochemically-coupled configurations. Based on these two simulations simplified expressions for the feedback parameters are obtained when the small temperature change in the biogeochemically-coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters allows identification of the reasons for differing responses among ocean and land carbon cycle models. </jats:p